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Abstract: The consequences for the brane cosmological evolution of energy ex- 
change between the brane and the bulk are analyzed in detail, in the context of 
a non-factorizable background geometry with vanishing effective cosmological con- 
stant on the brane. A rich variety of brane cosmologies is obtained, depending on 
the precise mechanism of energy transfer, the equation of state of brane-matter and 
the spatial topology. An accelerating era is generically a feature of our solutions. In 
the case of low-density flat universe more dark matter than in the conventional FRW 
picture is predicted. Spatially compact solutions are found to delay their re-collapse. 
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1. Introduction 

The idea that we might be living inside a defect, embedded in a higher dimensional 
space has already a long history. Concerning the nature of the defect, a solitonic 
codimension one or higher topological object was proposed [Q in the context of an 
ordinary higher dimensional gauge field theory, coupled || or not to gravity. It was 
soon realized, however, that, in contrast to scalar and spin-1/2 fields, it would be 
difficult to confine gauge fields on such an object. Various interesting ideas and 
scenaria were studied [[|, which have not yet given a fully satisfactory picture. In 
connection with the topology and the size of the bulk space on the other hand, 
the popular choice was that the extra dimensions are compact, with size of order 
O(Mpl). It was argued in ||, that in the context of the heterotic string with 
supersymmetry broken a la Scherk-Schwarz this should not be true anymore. The 
scale of supersymmetry breaking is tied to the size of internal dimensions, and a 
desirable supersymmetry breaking scale of a few TeV implies an extra dimension 
of about 10~ 16 cm. Despite difficulties to build a realistic model based on these 
ideas, the scenario was taken seriously and analyzed further for its phenomenological 
consequences ||. 
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The situation is drastically different in the context of type-I string theory ||. A 
few developments led to an exciting possibility and renewed interest in the whole idea. 
First, with the discovery of the D-branes as an essential part of the "spectrum" in 
type-I string theory, one could conjecture that we inhabit such a D-brane embedded 
in a ten-dimensional bulk. The usual solitonic defect of field theory was thus replaced 
by an appropriate collection of D-branes, which by construction confine the gauge 
fields 0], together with all the ingredients of the standard model. All known matter 
and forces lie on our brane world || ||, with the exception of gravity, which acts in 
the bulk space as well. It was, however, pointed out [K| , that for Kaluza-Klein extra 
dimensions the gravitational force on the brane was consistent with all laboratory and 
astrophysical experimental data, as long as the extra dimensions were smaller than 
a characteristic scale. This led to the exciting possibility of two extra dimensions 
in the sub-millimeter range. Furthermore, it was demonstrated in the context of 
an appropriate effective five- dimensional theory of gravity, that once we take into 
account the back reaction of the brane energy-momentum onto the geometry of 
space-time, the graviton is effectively confined on the brane and Newton's law is 
reproduced to an excellent accuracy at large distances, even with a non-compact 



extra dimension 11 



At the same time, the analysis of the cosmological consequences of the above 



hypotheses attracted considerable interest. The first step was taken in [12|, where 



the evolution of perfect fluid matter on the brane was studied, with no reference 
to the bulk dynamics and, therefore, no energy transfer between the brane and 
the bulk. Alternatively, a bulk-based point of view was adopted in where the 
cosmology induced on a moving 3-brane in a static Schwarzschild-AdSs background 
was studied, and a general interpretation of cosmology on moving branes, together 
with the idea of "mirage" cosmology, were presented in |14|]. The equivalence of 
the two approaches was demonstrated in [flj|. It was also pointed out that branes 



provide natural mechanisms for a varying speed of light [16 



Energy-exchange between the brane and the bulk should in principle be included 
in any realistic cosmological scenario, and its effects have been studied in detail in 
the context of flat compact extra dimensions [17]. The role of energy-exchange on 
brane cosmology in the case of non-factorizable extra dimensions has not yet been 
investigated extensively [Tj| |TU[, even though the importance of energy outflow from 
the brane has already been demonstrated in the context of a Randall- Sundrum 
configuration with additional gravity induced on the brane . 

The present paper is an attempt towards a more complete analysis of the cosmic 
evolution of the brane in the presence of energy flow into or from the bulk. Our aim 
is to generalize the picture in bulk AdS (bulk gravity plus cosmological constant) by 
considering a general bulk theory (that includes gravity). There may be more bulk 
fields and more general bulk-brane couplings. Although this can be formulated by 
the standard action principle and the relevant exact equations derived and studied 
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we have opted in this paper for a short-cut. We are analyzing the regime where the 
bulk energy, at the brane position, can be consistently neglected from the equations. 
Moreover, we parameterize appropriately the energy exchange as a specific power of 
the matter density of the brane. Although there are other valid parameterizations, 
we opted for this one, motivated by our previous work in []nj where we analyzed the 
out-flow of energy from the brane due to graviton radiation in the presence of induced 
gravity. Standard cross section calculations give an out-flow that is a function of the 
temperature and other fundamental constants of the theory. If one re-expresses the 
temperature in terms of the running density using the cosmological equations we 
obtain a rate of flow that is a power of the density, with a dimension-full coefficient 
that depends on fundamental constants and initial energy densities. This is valid 
typically for a whole era, that is, piecewise in the cosmological evolution. 

In a given theory, with a specific bulk content and brane-bulk couplings, the 
exponent of the density as well as the coefficient are calculable functions of the 
coupling constants of the Lagrangian as well as initial densities. It is also obvious 
that if the bulk theory is approximately conformal the most general form of rate of 
outflow will be polynomial in the density. Our analysis is "phenomenological" and 
we parameterize different theories with different forms for the energy in/out-flow. 

The paper is organized in five sections of which this introduction is the first. 
In section 2 the framework of our work is described and the approximations on the 
brane-bulk exchange are presented. The effective equations for the analysis of the 
brane cosmology are derived. Section 3 contains several interesting characteristic 
solutions of the brane cosmology, while an exact treatment of the influx/outflow 
equations, relevant for a wide range of potentially realistic applications, is presented 
in section 4. Our results are summarized and the prospects for further research along 
these lines are discussed in the final section. 



2. The model 

We shall be interested in the model described by the action 

S = Jd 5 x ^=g~ (M Z R - A + C™ at ) + J d 4 x^g~ (-V + , (2- 1 ) 

where R is the curvature scalar of the five- dimensional metric Qab, A,B — 0, 1, 2, 3, 5, 
A is the bulk cosmological constant, and g a/ 3, with a,/3 = 0,1,2,3, is the induced 
metric on the 3-brane. We identify (x, z) with (x, — z), where z = x§. However, 
following the conventions of |]TTJ we extend the bulk integration over the entire in- 



terval (— oo,oo). The quantity V includes the brane tension as well as quantum 
contributions to the four- dimensional cosmological constant. 
We consider an ansatz for the metric of the form 

ds 2 = -n 2 (t, z)dt 2 + a 2 (t, z)j ij dx i dx j + b 2 (t, z)dz 2 , (2.2) 
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where 7^ is a maximally symmetric 3-dimensional metric. We use k = —1,0,1 to 
parameterize the spatial curvature. 

The non-zero components of the five- dimensional Einstein tensor are 
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Primes indicate derivatives with respect to z, while dots derivatives with respect to 
t. 

The five- dimensional Einstein equations take the usual form 

1 
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(2.7) 



where Tac denotes the total energy-momentum tensor. 

Assuming a perfect fluid on the brane and, possibly an additional energy-momentum 
Tc\m,B in the bulk, we write 
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diag(-\/,-F,-K,-^,0) 



diag(-A, -A, -A, -A, -A) 
—r- diag(-p,p,p,p, 0), 



(2.8) 

(2.9) 
(2.10) 

(2.11) 



where p and p are the energy density and pressure on the brane, respectively. The 
behavior of T^\ m ^B is in general complicated in the presence of flows, but we do not 
have to specify it further at this point. 

We wish to solve the Einstein equations at the location of the brane. We indicate 
by the subscript o the value of various quantities on the brane. Integrating equations 



D' (HH) with respect to z around z = gives the known jump conditions 
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The other two Einstein equations (|2.5| ), (|2.6|) give 

a 'o bo = 1 

a D ~ 6M 3 
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(2.15) 



where T 05 ,T 55 are the 05 and 55 components of TAc\m,B evaluated on the brane. 
Substituting ( gTp , (|2~l"3l) in equations (gH , (gl5|) one obtains 



p + 3— (p + p) 
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(2.16) 





(A 



1 



144M 6 



6M 3 V 12M 3 
(V(3p-p)+p(3p + p)) 



T"l5 

6M 3 5 ' 



(2.17) 



We are interested in a model that reduces to the Randall- Sundrum vacuum |[T | 
in the absence of matter. In this case, the first term on the right hand side of equation 
( |2.17D vanishes. A new scale k^s is defined through the relations V = —A/k^s — 
12M s k RS . 

In order to derive a solution that is largely independent of the bulk dynamics, the 
T\ term on the right hand side of the same equation must be negligible relative to the 
second one. This is possible if we assume that the diagonal elements of the various 
contributions to the energy-momentum tensor satisfy the schematic inequality 1 



rp I diag 
1 lm,B 



rp I diag 
1 \y,B 



rp I diag 
lm,6 



rp I diag 
1 \v,b 



(2.18) 



Our assumption is that the bulk matter relative to the bulk vacuum energy is much 
less important than the brane matter relative to the brane vacuum energy. In this 
case the bulk is largely unperturbed by the exchange of energy with the brane. 
When the off-diagonal term T° 5 is of the same order of magnitude or smaller than 
the diagonal ones, the inequality ( [2.18| ) implies T pk^s- 

At this point we find it convenient to employ a coordinate frame in which 
b — 7i o = 1 in the above equations. This can be achieved by using Gauss nor- 
mal coordinates with b(t,z) = 1, and by going to the temporal gauge on the brane 



1 Strictly speaking, the left hand side of (|2.18| ) concerns only the 55 components of the bulk 
contributions to the energy-momentum tensor. The other components do not appear in equations 
( 2.16| ), (2.17) and do not affect the cosmological evolution on the brane. 



5 



with n — 1. The assumptions for the form of the energy-momentum tensor are then 
specific to this frame 2 

Using (3 = M _6 /144 and 7 = V r M _6 /144, and omitting the subscript o for 
convenience in the following, we rewrite equations ( |2.16| ) and (|2.17 ) in the equivalent 
form 

p + 3(1 + w) - p = —T (2.20) 
a 

- 2 =f3p 2 + 2 1P -^ + X + \ (2.21) 
a 2 a 2 

(2.22) 

where p = wp, T = 2T° 5 is the discontinuity of the zero- five component of the bulk 
energy-momentum tensor, and A = (A + V 2 /12M 3 )/12M 3 the effective cosmological 
constant on the brane. 

In the equations above, Eq. (|2.21|) is the definition of the auxiliary density x- 
With this definition the other two equations are equivalent to the original system 




fl2.16|j2T7|) . As we will see later on, in the special case of no-exchange (T = 0) x 
represents the mirage radiation reflecting the non-zero Weyl tensor of the bulk. 
The second order equation (p,17|) for the scale factor becomes 



- = -(2 + 3w)pp 2 - (1 + 3w) 1P - x + A. (2.23) 
a 

As mentioned above, in the Randall- Sundrum model the effective cosmological con- 
stant A vanishes, and this is the value we shall assume in the rest of the paper. 

In the special case of w = 1/3 one may define a new function x = X + ^1P- The 
functions x, p and a satisfy equations (|2.20| ) to ( |2.23| ) with x i n place of x an d 7 = 0. 



This should be expected, since for w = 1/3 there is no 7 left in equation ( [2.17|) . 



3. Special solutions 



Before presenting exact solutions of equations ( |2.20| )-( ]2.22| ), it is instructive to con- 
sider a few special cases whose physical content is more transparent. 



2 If the vacuum energy dominates over the matter content of the bulk, we expect that the form 
of the metric will be close to the Randall-Sundrum solution with a static bulk. Thus, we expect 
(even though we cannot demonstrate explicitly without a full solution in the bulk) that in a generic 
frame, in which 
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we shall have b ~ 0. Then the transformation that sets b - 
the energy- momentum tensor. 



(2-19) 



1 is not expected to modify significantly 
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We concentrate on the low-density region, in which p <C 7//?. In this case we may 
ignore the term (3p 2 in the above equations compared to jp. As a result, equations 
(gggHggp can be written as 



p + 3(1 + w)Hp = — T 



-)' 

ay 



k 

27P + X - -o 
a 2 



(3.1) 

(3.2) 
(3.3) 



The cosmological evolution is determined by three initial parameters (pi, a i: Xh or 
alternatively p^ hi), instead of the two (pj, cjj) in conventional cosmology. The 



reason is that the generalized Friedmann equation (|2.21|) (or |3.2| ) is not a first integral 
of the Einstein equations because of the possible energy exchange between the brane 
and the bulk. 



3.1 "Mirage" radiation for energy outflow 

Let us consider first the case T > 0, for which there is flow of energy out of the brane. 
If the brane matter is radiation dominated with p = p/3, equations fl3.1|), ( p.3[ ) have 
an exact solution independently of the explicit form of T: 



P + 



X_ 

2 7 



Pi + 



27, 



(3.4) 



and 



H* = (2 1 p t + Xl )- A -- 2 . (3.5) 
a 4 a 1 

Assume that initially Xi = 0- It is clear that the effect of the radiation on the 
expansion does not disappear even if it decays during the cosmological evolution: 
the Hubble parameter of equation (|3.5| ) is determined by the initial value of the 
energy density, diluted by the expansion in a radiation dominated universe. The real 
radiation energy density p, however, falls with time faster than a -4 . 

As a simple example we consider T = Ap with A > 0. Then equation (|3.1|) can 
be integrated, with the result 



P = Pi 

where we have considered the general case p 
Xi = 0, w = 1/3 we obtain 



X 



-At 



wp 
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(3.6) 



From equation (3.4) with 



(3.7) 



3 Having neglected the term j3p 2 from equations (3^) and (B^), these solutions are valid only 
if (3p1 <C Hf. As we shall show in subsection (4.2) this condition is eventually satisfied for all 
solutions in the case of outflow of the form discussed here and k = 0. 
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The Hubble parameter, given by equation with Xt — 0, corresponds to an 

initial radiation density pi, further diluted only by the expansion. At late times the 
expansion is the consequence of a "mirage" effect, arising from the original radiation 
through x- 

The presence of a "mirage" term is possible even without energy flow. More 
specifically, equation ( |3.3|) has a solution x — C /a A even for T = 0, and x can ac t as 
"mirage" radiation. The novel feature for T 7^ is that the "mirage" effect appears 
through the decay of real brane matter, even if it was absent in the beginning. 

The effect persists even if the decaying matter is not radiation. Consider a 
more realistic case, in which the brane matter consists of non-relativistic particles 
with p = 0, that can decay spontaneously into bulk matter. In this case our previous 
assumption T = Ap becomes realistic. The brane energy density is given by equation 
( |3.6| ) with w = 0. For Xi — we find 



For t — > 00, if the scale factor grows more slowly than an exponential, the integral 
converges and % ~ 1/a 4 . Again, there is a "mirage" effect corresponding to an initial 
energy density proportional to pj. 

It can be checked easily that, if the decaying matter is only one component of the 
matter on the brane, the "mirage" effect acts as an additional radiation component 
at late times. Whether it can affect the density perturbations of real matter and act 
as hot dark matter is an open question. 

For the solution (|3.4j), the quantity q = a/a equals —{2^pi + Xj)a 4 /a 4 = —a/a 4 . 
For k > 0, equation (|3.5|) requires a > and, therefore, q < 0. For k = — 1, o may 
be negative, leading to a loitering universe. The radiation case is not realistic for 
the eternal stage of the universe, but we are going to see in the following a situation 
with k = — 1, w = and eternal acceleration. 

3.2 Accelerating solutions 

An interesting feature of this framework is the possible presence of accelerating cos- 
mological solutions. We can look for exponential expansion with a constant Hubble 
parameter H, even if the brane content is not pure vacuum energy. This implies 
that equations ( |3.1|) - (|373|) must have a time-independent solution, without necessar- 
ily requiring w = — 1. The possible fixed points (denoted by *) of these equations 
for k = satisfy 






(3.9) 
(3.10) 
(3.11) 



Hi = 27P* + x* 
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It is clear from equation (|3T9|) that, for positive matter density on the brane (p > 0), 
flow of energy into the brane (T(p) < 0) is necessary. 

The accretion of energy from the bulk depends on the dynamical mechanism that 
localizes particles on the brane. Its details are outside the scope of our discussion. 
However, it is not difficult to imagine scenaria that would lead to accretion. If the 
brane initially has very low energy density, energy can by transferred onto it by bulk 
particles such as gravitons. An equilibrium is expected to set in if the brane energy 
density reaches a limiting value. As a result, a physically motivated behavior for the 
function T(p) is to be negative for small p and cross zero towards positive values for 
larger densities. In the case of accretion it is also natural to expect that the energy 
transfer approaches a negative constant value for p — > 0. 

The solution of equations (|3.9|) - (|3.11|) satisfies 



T(p,) = - 3 -^-(l + w)(l-3w)^p 3 J 2 (3.12) 
Hi = 1 -^1P* (3.13) 
X, = (3.14) 

For a general form of T(p) equation ( |3.12| ) is an algebraic equation with a discrete 
number of roots. For any value of w in the region — 1 < to< l/3a solution is possible. 
The corresponding cosmological model has a scale factor that grows exponentially 
with time. The energy density on the brane remains constant due to the energy flow 
from the bulk. Our model is very similar to the steady state model of cosmology [f2~0| 



The main differences are that the energy density is not spontaneously generated, and 
the Hubble parameter receives an additional contribution from the "mirage" field \ 
(see equation ( |3.10| )). 

The stability of the fixed point ( |3.9| )- (|3~Tl| ) determines whether the exponentially 
expanding solution is an attractor of neighboring cosmological flows. If we consider 
a small homogeneous perturbation around the fixed point (p = p* + Sp, x = X* + $x) 
we find that Sp, 8\ satisfy 



(3.15) 



where 

/- v + 3(1 - w)/(l - 3w) 1/7(1 - 3w) \ 

{ 2 T (i? - 2/(1 - 3u>)) -2(l + 9w)/[3(l + w)(l-3w))J 1 ! 





d\n\T\ 
dlnp 



(p*) , (3-17) 
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and we have employed the relations (|3.9| )- (|3.14|) . If the energy flow has a simple 
form T(p) oc p u , we have v = v. The eigenvalues of the matrix M. are 



_ 7 + 3w - 3z>(l + w)± V24(-3 + 2x/)(l + iu) + [7 + 3w - 3z>(l + w)] 2 

6(1 + w) 

(3.18) 

For — 1 < w < 1/3, < v < 3/2 they both have a positive real part. As we have 
assumed T(p) < 0, the fixed point is stable in this case. We have verified this result 
by integrating equations (|3.1|)-Q3~3|) numerically. The approach to the fixed-point 
values depends on the sign of the quantity under the square root. If this is negative 
the energy density oscillates with diminishing amplitude around its fixed-point value. 

For w = — 1 we get the standard inflation only for a value p* that is a zero of 
T(p). In this case there is no flow along the fifth dimension and also x* = 0- 

3.3 The case of radiation for energy influx 

In the case of radiation the general solution of equations (|3.1|)-(pT3|) was derived 
in the previous subsection and is given by equations (|3.4j ), (|3~5|). The expansion is 
that of a radiation-dominated universe with constant energy (xi/{2j) + Pi)a\ per 
co-moving volume. The "mirage" energy density is diluted ~ a -4 . 

The explicit dependence on time will be discussed next in the case of flat space 
(k = 0), in which the energy density satisfies 

| + f, = -T(,). (3,9) 

If T(p) < for all p, and the "friction" term in the left hand side becomes suppressed 
for t — > oo, we expect an unbounded increase of p in this limit. For p ^ j//3 the 
low energy approximation employed in this section breaks down. The full treatment 
necessary in this case will be given in the next section. 

The actual situation is rather complicated and the details depend crucially on 
the form of T(p). Assuming that T(p) = Ap v with A < 0, the exact solution of 
equation ( |3.19| ) for v ^ 1, 3/2 is 

pp) lu = ( P m Xv + 1^- ft 3 -" - r , (3.20) 



2v 

where t = \A\t. For v — 1 the solution is 



and for v = 3/2 



t 2 - - 

P = Pij-J~ t \ (3-21) 




(3.22) 
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For < v < 1 we have p ~ for t — > oo. For v — 1 the increase of the 

energy density at large t is exponential moderated by a power. For 1 < v < 3/2 the 
energy density diverges at a finite time 

^ = ^ + ^(^1)^ (3-23) 



1-u 

A similar divergence appears for v = 3/2. For v > 3/2 a divergence occurs if the 
quantity 

^ = (p^)^-^^- 3 (3.24) 

is negative. In the opposite case pt 2 — > 1/D for t — > oo, and the energy density 
diminishes: p ~ t~ 2 . 

As we discussed earlier, it is physically reasonable that the energy influx should 
stop at a certain value p cr , and be reversed for larger energy densities. The dynam- 
ical mechanism that localizes particles on the brane cannot operate for arbitrarily 
large energy densities. This modifies the solutions above that predict an unbounded 
increase of the energy density. 

A final observation that will be encountered again in the next section is that, 
despite of the fact that the energy density in most cases increases for large times, it 
can decrease at the initial stages. This is obvious from eq. ( [3.19|) . If at the time t r 
that the brane enters a radiation dominated era \T(p)\/p < 2/t, the energy density 
decreases for a certain time. 

3.4 Non-flat solutions 

In addition to the analytical special solutions discussed above, we would like to 
present a few suggestive numerical results concerning the k = ±1 cases. For v — 1, 



we substitute ( |3.6|) , true for any A, into the second order equation fl2.17|) , to obtain 



I + £ + £ + (i + MM - (i - *»h* = o • < 3 - 25 > 



It is obvious from ( |3.25|) that for outflow with k — 1 and w > —1/3 the Universe will 



exhibit eternal deceleration. In particular, in the case of dust, Figure 1 depicts the 
solution for a(t) of ( |3.25| ) for some initial conditions for a, p. Notice that a(t) after a 



period of decrease, starts increasing again, thus delaying its re-collapse. Of course, 
with appropriate initial conditions one also obtains solutions with the standard FRW 
behavior. 

For an open (k = —1) Universe with A > 0, w = 0, v = 1 we have found 
numerically a solution where p(t) monotonically decreases to zero, while a(t) starts 
with deceleration, but later on accelerates eternally; more specifically, q(t) — > + for 
t — > 0. Another possibility for k = —1 and A < 0, w = 0, u = 1 allows for a Universe 
starting with acceleration at infinite densities, which later on, turns to deceleration 
with p approaching a constant value. 
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a 







Figure 1: Outflow, k = +1, w = 0, v = 1. 

arguments: 

a) The evolution equations for T = 



4. General treatment of the out- 
flow/influx equations 

In this section, we present a few first steps 
towards an exact analytical treatment of 
equations (^20|) -(^22D that govern the cos- 
mological evolution. As we mentioned ear- 
lier, we concentrate on the case with zero 
effective cosmological constant A and in- 
vestigate the brane cosmological dynamics 
due to the energy exchange with the bulk. 

Some qualitative properties of the ex- 
act solutions can be derived from general 



reduce to those of ref. jl2| : ( |2.20|) is the 
standard energy conservation equation, while fl2.22j ) upon integration leads to the 
usual mirage radiation term \ = C/a A . If there is no energy exchange between the 
brane and the bulk (T = 0), it can be seen from equation ( |2.23|) that for suitable 
C < an acceleration era (a > 0) is possible. It is tempting to associate this era 
with the recently observed cosmic acceleration JH]. It is an open question, however, 
if the preceding deceleration era (with a dominant contribution ~ p 2 to the Hubble 
parameter) can accommodate the successes of standard Big-Bang cosmology. 

b) For T^O, a suitable negative x term can, in principle, produce acceleration 
(i.e. behave like a positive cosmological constant in equation fl2.23| )) and at the same 
time play a role similar to that of a negative cosmological constant in the Hubble 
expansion ( |2.21|) . A consequence of this fact is that, if equations Q2.20|) -( f2.23| ) are 
to describe the present-day Universe with a/ a > 0, we must require x < 0. In 
the asymptotic regime p <C t//3, this means that, for k = 0, the matter content of 
the critical Universe must compensate for the negative contribution from \. This 
requires fl m > 2 and Q x < — 1. As a result, this scenario predicts more dark energy 
than conventional FRWL models. Notice, however, that for k = — 1 these restrictions 
do not apply. For example, we could have Qk = 0.96, Q x = —0.06 and fl m = 0.1. 

c) In the case of outflow (T > 0), it is obvious from the conservation equation 
( |2.20p that p(t) is monotonically decreasing for all expanding solutions. Furthermore, 
if the function T(p) vanishes only at p = 0, p(t) — > for t —* +oo. 

d) It is straightforward to argue that for k = 0, +1 and w > —1/3, when p > 
(1 — 3w)j/((l + 3w)P) we have necessarily deceleration. For k = 0, +1, w > 1/3 
we have deceleration at all times. The absolute upper bound for a /a is at all times 
(1 — 3w) 2 7 2 / (4(1 + 3w)/3). These conclusions follow from the fact that the left hand 
side of equation ( |2.2 1| ) must be positive definite. This implies that a/ a lies below the 
parabola (1 — 3w)jp — (1 + 3w)j3p 2 , which for w > 1/3 is negative definite. Finally, 
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one may see from equations ( 2.21|) , (|2.22 ) that all expanding k = solutions with 
p — > for t — > +00, will have x ~ > 0, d/a — > CT and d/a — > + , for functions T(p) 
with T(0) = 0, while for T(0) < 0, p — > is not possible. 

4.1 General analysis of the equations 

Combining equations ( p. 20 ) and ( 2.21 ) we obtain 

dp 



-1/2 



da 



-3(1 + w)p - eT(p) /V + 2 7 p 



+ X 



(4.1) 



Similarly, equations (|2.22j) and (|2.21|) give 



27P - + X 



'1/2 



(4.2) 



where e = 1 refers to expansion, while e = — 1 to contraction. These two equations 
form a two-dimensional dynamical system. The function x{p) is obtained from the 
equation 



k 



dx 



3(1 + w)p^/3p 2 + 2 7 p - - + X + eT(p) 1 -f 



2e/5 (p + ^j T(p). 



(4.3) 



Note that the equations of contraction are those of expansion with the roles of out- 
flow and influx interchanged. Combining equations ( |4.3j ) with ( |2.23| ) we derive the 
following equation for q(p): 

dq{p) 

dp 



3(1 + w)p + eT(p)/J(l - Zw) 1P - (1 + 3w)(3p 2 - ka~ 2 - q(p) 



+e(2(l + 3w)(3p - (1 - 3w) 1 )T(p)/^J(l - 3w) 1P - (1 + 3w)f3p 2 - ka~ 2 - q(p) 
Aq{p) + (1 + 3w) (2(2 + 3w)pp - (1 - 3wh)p = . (4.4) 



We shall use this equation below to decide about eras of acceleration and deceleration 
in the cosmic evolution on the brane. For \w\ < 1/3, T(p) > and p < ((1 — 
3w)7)/(2(2 + 3w)/3), it can be seen from the above equation that q > implies 
dq/dp > 0. However, in the range ((1 - 3w)7)/(2(l + 3w)(3) < p < ((1 -3w) 7 )/((l + 
3w)P), passage from acceleration to deceleration during expansion is not allowed. 

Another form of the system of equations ( |4.1| ), ( |4.2j ), which will be helpful for 
the study of fixed points of the system, can be derived by defining a = In a and 
Q = q + ke~ 2a . Then 



da 



2a 



+ S(p,Q), 



(4.5) 
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where 

R(p, Q) = -3(1 + w)p - eT[-{\ + 3w)(3p 2 + (1 - 3w) 1P - Qy l/2 } (4.6) 
S(p, Q) = (1 + 3u>)(2(2 + 3w)/3p - (1 - 3w) 1 )p - 4Q 

+eT(2(l + 3w)f3p - (1 - 3w) 7 ) (-(1 + 3w)pp 2 + (1 - 3w) 7 p - Q) ~ V2 . 

(4.7) 

When the 3-metric is flat (k = 0), equation (|4.3|) becomes an autonomous equa- 
tion, which, in principle, may be solved for x(p)- We define 

2 / d 5+3w 

Y(p)=p sa-H.) ^/5p2 + 2 7p + x + —— -p W+^T(p), (4.8) 

v 3(1 + w) 

in terms of which equation (|4.3|) translates to 



dY e 8+6^ / 7 + 3w m dT\^ r 2 13+9^ 

K — -I n m+w) [ T — p — \Y p w+^jt 

dp 3(l + w) \3{l + w) H dpJ 27(1 + «; Na 



2 



3' 



{ l-3w 1 + 3w \ 4 / . 

+ \3(T+wj' T ~ 3(1+10)^7 P = (49) 

We shall simplify our discussion by considering energy transfers of the form 

eT( P ) = Ap\ (4.10) 

with 1/ an arbitrary real parameter. The brane-bulk energy exchange in the real 
world may be much more complicated. For example, it may be a sum of terms of 
the above form. 

With the change of variable 

p = p u W+^, (4.11) 

the equation for Y becomes 

Y d 4- - BY 
dp 

1 / 2A 2 \ 

5 + 3w — 3v{l + w) y 9(1 + w)^ / 

(4.12) 

where 5 = A(7 + 3w - 3z/(l + w))/(3(l + w){5 + 3w - 3u(l + w))), r = (7 + 9w - 
3u{l + w))/(5 + 3w-3u(l+w)) and s = {A + 6w - 3u(l +w))/{5 + 3w - 3u{l + w)). 

In the cases with — 1 < w < 1/3, v < (4 + 6u>)/(3(l + u>)), or w > 1/3, 
1/ < (5 + 3w)/(3(l + w)), or -1 < w < -1/3, 1/ > (5 + 3tu)/(3(l + tu)), or w > -1/3, 
1/ > (7+9u/j/(3(l+w)), we haver, s > 0. In the case w > 1/3, (5+3w)/(3(l+tt;)) < 
v < 3/2, we have r, s < —1. In all these cases, the behavior p « 0, if such a region 
exists, is given by the following approximation of equation ( [4.12 ): 



dY 2A 2 

Y—-BY + - — ttP = 0. (4.13) 

dp 9(l + w) 2 (5 + 3w -3v(l + w)) H V ; 
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4.2 Solutions with p « 

Case i: For v = \ (and thus to > —1/3), as in the case of unstable matter on the 
brane decaying into the bulk, the general solution of (|4.13|) is |22 



(-^M^h^ (4 - 14) 

Equivalently, in terms of the function x, the above solution takes the form 

A 1 \ a 



(3p 2 + 2 7 p + x exp - = k P WT^ } (4.15) 

V 3(1 + w) y /Pp i + 2n/p + xJ 

where k is a non- negative integration constant. Using (|2.21[ ) we may rewrite the 
above as 

|#| e ~3(iT^WT = Kp3^+™>. (4.16) 

Equation (|4.16|) , as is suggested by its independence on (3 and 7, is valid only 
as long as the (3 and 7 terms in ( |4.12|) are negligible compared to the linear term 
inside the parenthesis. In the case of outflow this is not a constraint, since p is 
monotonically decreasing towards zero. So, at late times all solutions are correctly 
described by ( f4.16|) . 



From equation ( |4.14j ) we conclude that for A < (influx), v = 1, w > — 1/3, any 



solution p(t) with k > is bounded from below by some positive value. This implies 
that the point (p = 0, q = 0) cannot be an attractor. Furthermore, using Q4.16| ) and 
Q2.17D one may show that for w ^ 1/3 q ^ (1 — 3w)'j P , i.e. all curves near the origin 
of the (q, p) plane have a fixed slope. This will be confirmed numerically in Figure 3 
below. Also, \H\ < p/M 3 . 

In the case of outflow, on the other hand, equation ( 4.16|) implies that for v = 



1, w > —1/3, we have \H\ — > for p — ► 0. Equation ( |4.16| ) implies that p \q\ and 



in addition that p/q ^ — ooasp^O, in agreement with the case depicted in Figure 
2. Using ( |2.20|) - ([2.22|) it is straightforward to show that 

dH 2 

— — + AH 2 - 2(1 - 3w) 1P + 2(1 + 3w)fip 2 = 0, (4.17) 
da 

which combined with (|4.16| ) leads to 



\H\ = ~ 2 , (4-18) 
cr 

or equivalently to 

a 2 {t) = 2 Cl t + c 2 , (4.19) 

with ci and c 2 integration constants. Using this expression of H(t) in ( |4.16|) , we 
obtain 

p(t) = (2.)- 3 ^ 2 (t+ e ^ )3(1+w)/2 . (4.20) 
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Apart from the obvious agreement of the above with the special solution (|3.21| ) (also 
valid for A > 0), we should like to stress at this point the remarkable fact that the 
scale factor a(t) of the Brane Universe, for v = 1 and essentially any w, behaves at 
late times as if radiation dominated. 

A global phase portrait of q = a /a 
with respect to p during expansion in the 
outflow case for k = 0, w = 0, v = 1 
is shown in Figure 2. The characteris- 
tics of the solutions are in agreement with 
our conclusions above. One recognizes two 
families of solutions: The first have q < 
for all values of p, while the second start 
with a deceleration era for large p, enter 
an acceleration era and then return to de- 
celeration for small enough values of p. 
This is exactly what should be expected 
on the basis of the previous analysis. So- 
lutions corresponding to initial conditions 
with positive q (always under the limiting parabola shown with the dotted line), 
necessarily had a deceleration era in the past, and are going to end with an eter- 
nal deceleration era also. The straight dashed line represents the standard FRW 
solution. 

Case 2: For v ^ 1, following the same steps as in case 1, the general solution of 
(gT3) for -l<w< 1/3, v < (4+6w)/(3(l+w)), orw> 1/3, v < (5+3w)/(3(l+w)) 




Figure 2: Outflow, k = 0, w = 0, v = 1. 

The arrows show the direction of increasing 
scale factor. 



is 



Ki 



Y 



A 



3(l + w 



-P 



5+3to-3i/(1+id) 
2 



«2 



Y - 



2A 



3(1 + w)(5 + 3w - 3v(l + w)) 



P 



(4.21) 



where |ki| + | At 2 1 > 0. 
solution for H(p) 



Using ( |4.8| ) and ( |2.21[ ) this is translated into the implicit 



/Cilifl 



5+3to-3i<(1+iu) 
1 



«2 



H\ - 



A{v-l) 



5 + 3w - 3z/(l + w) 



(4.22) 



A few comments concerning the validity of our approximations are in order at 
this point. First, we would like to stress that ( |4.14|) and fl4.21|) were derived and are 
valid only near "vanishing" p. Their validity may be questionable in situations with p 
bounded from below. A detailed perturbative analysis of the next order correction to 
Y(p), ensures that the above expressions apply in all cases studied, except for some 
sets of initial data in the case of energy influx domination (A < 0, v < 1), in which 
p is indeed bounded from below. A similar comment applies to equations ( |4 . 1 6| ) and 
( f4.22| ). The translation from Y to H involves equation ( |4.8| ), whose possible singular 
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nature leads to some further constraints on their validity. Specifically, in the case 
of v > 1 for both influx and outflow, ( |4.16|) and ( |4.22j) are correct for all solutions 
with p approaching arbitrarily close to zero. For dominant outflow or influx, i.e. 
v < 1, ( f4.22| ) should not be trusted. In these cases, the next order correction to H 



is of the same order of magnitude as the leading if, given in ( |4 . 2 2| ) . Incidentally, 
the above clarifications remove an apparent contradiction of ( |4.22j ) with the general 
comment (d) in the beginning of this section. Namely, ([4.22|) implies that for v < 1, 
H(p) tends to a non-vanishing constant as p —>■ 0. This contradicts comment (d) but 
happens in a region of parameters and initial data for which ( |4.22|) is not reliable. 

Finally, for 1 < v < 3/2 it is straightforward to see that q/ p — > — oo and, 
furthermore, in the case of outflow, the explicit expression for the scale factor is 
given by ( |4.19| ), while the matter density is 

/ c 1 (5 + 3w-3u(l+w))\ 1/{u ' 1) 



Pit) V A(u-l)(2 Cl t + c 2 ) J 
In the case of influx the point (0, 0) in the (p — q) plane is a repulsor. 



(4.23) 



4.3 Fixed points 

Equations ( |2.16| ) and ( j2.17| ) possess for k = the obvious fixed point solution (p* = 
0, if* = 0). However, there are more fixed points which may be found by setting 
dp/ da = dQ/da 
relations 



in the system of equations fl4.5|) . This leads to the following 

2|T(p*)| 2 



3w)Ppl 



q* 



3w)7p* 
\T(P*)\ 2 



9{l + wYpl 



> 0. 



0. 



(4.24) 



(4.25) 



9(l + u>) 2 p^ 

It is obvious from the conservation equation ( |4.1|) that during expansion non-trivial 
fixed points may exist only in the influx case. Their number is determined by the 
roots of equation ( 4.24 ). It can be seen from ( 4.25| ) that these solutions are acceler- 
ating. From ( p^ID , (ggg ) we find that = ^ and X * = -3(1 + w)(@p* + 7)p*/2. 

In order to study the stability of these fixed points we write p(a) = p* + Sp(a), 
q(a) = g* + 5q(a). The resulting linearized equations are 

d ( 5p(a) \ _ 9(1 + w) 2 pl { m 1 m 2 \ ( 5p(a) 
da \ 5q(a) J 

where 



2T 2 



mi m 2 

TO3 m.4 



5q{a) 



(4.26) 



mi = 3(1 + w) ((3 - z/)(l + 3w)p Pif + {v - 2)(1 - 3w)j) 
m 2 = 3(1 + w) 

m 3 = ((1 - 3w)-f - 2(1 + 3w)pp,) 

x [(1 + 3w){7 + 9w- 3z/(l + w))0p* - (1 - 3w){4 + Qw - 3u{l + w))j] 



m 4 = -(2(1 + 3w) 2 (3p* + (1 - 3w) 2 7 ). 



(4.27) 
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The sign of the real part of the eigenvalues of the matrix appearing above determines 
the nature of the fixed point. The presence of an imaginary part results in a spiral 
form for the flows. 

As an example, we consider the case with dust (w = 0) and v — 1. For \A\ < 
3^/^/8/3, there exist two real positive roots 

The root pl +) is always a saddle point. Furthermore, for 27 / y 2 /(32/3) < A 2 < 
97 2 /(8/9), pi~^ is a stable node, while for A 2 < 277 2 /(32/3), it becomes a coun- 
terclockwise stable spiral. For (3 A 2 /^ 2 < 1 we have pi ■* ~ 2A 2 /(9'y) < 7//? and 
jO* ~ lift- In this limit, pi ^ corresponds to the fixed point ( |3.12| ) of the spe- 
cial solutions discussed in section 3. The eigenvalues of the stability matrix are in 
agreement with the ones derived there. 

The global phase portrait of q = a/ a with respect to p during expansion for 
the case k = 0, w = 0, v = 1 is shown in Figure 3. The presence of the limiting 
parabola as in the outflow case is apparent. However, new characteristics appear. 
For example, pi - "* attracts to eternal acceleration a whole family of solutions which 
start their evolution at either very low or very high densities. There is another family 
of solutions which are attracted to acceleration by p{ + ^ and which eventually exit to 
a deceleration era. Finally, there is a family of solutions, near the limiting parabola, 
which start with acceleration at very low densities, and eventually exit to eternal 
deceleration, while their density increases monotonically with time because of the 
influx. 




Clearly, for v 7^ 1 one expects a dif- 
ferent set of fixed points with varying be- 
haviors around them. 

5. Conclusions 

A rather detailed mathematical analysis 
of the role of the brane-bulk energy ex- 
change on the evolution of a Brane Uni- 
verse, adequate for a wide range of poten- 



Figure 3: Influx, k = 0, w = 0, v = 1. tiall y realistic implementations, was pre- 

sented. The effective brane cosmological 
equations were derived with perfect fluid matter on the brane, constant energy- 
momentum tensor in the bulk and non-vanishing exchange between them. A de- 
tailed study of the solutions of these unconventional equations was performed in the 
case of zero effective cosmological constant on the brane, in order to reduce to the 
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Randall- Sundrum vacuum in the absence of matter. The analysis revealed a rather 
rich variety of possible cosmologies, depending upon the precise form of the exchange 
term, the topology of 3-space, and the nature of matter on the brane. 

A few special but particularly interesting solutions were obtained in the limit 
of low energy density on the brane. One of these is the exactly solvable case of 
radiation, where a mirage radiation effect appears through the decay of real brane 
matter. Another, is a De Sitter fixed point solution, obtained in the case of energy 
influx, even without pure vacuum energy, which in addition is stable for a wide range 
of reasonable forms of energy exchange. 

An exact treatment of the full set of equations followed. Several generic qual- 
itative bounds were extracted and the asymptotic characteristics of the solutions 
were studied. In particular, in the low-density regime and for a flat universe, our 
model predicts additional dark matter compared to the standard cosmological pic- 
ture. Adopting the physically motivated T ~ p u power-law form for the energy 
transfer, we found exact analytical solutions in the low-density regime for almost all 
values of v and of the matter equation of state. Furthermore, the global phase por- 
traits in the density-acceleration plane were obtained numerically, separately for the 
outflow and the influx cases, for arbitrary values of the density. In the outflow case, 
one distinguishes two families of solutions: one describing an all-time decelerating 
universe, and another which describes a Universe with an intermediate accelerating 
era. The latter starts decelerating, enters an accelerating phase, and finally deceler- 
ates again, with the energy density decaying to zero. In the case of influx on the other 
hand, one generically obtains several fixed point solutions, whose stability analysis 
was performed. Curiously, all these fixed points correspond to positive acceleration, 
and their presence implies the existence of solutions with accelerating era. In addi- 
tion, as one may see in Figure 3, there are solutions that can be purely accelerating, 
purely decelerating, or with alternating accelerating and decelerating behavior. 

A complete description of the evolution of our Brane Universe is still lacking. 
Several scenaria are possible and their detailed study in the light of observations will 
determine their viability. For example, a sketchy cosmological evolution could be as 
follows: The brane is created much hotter than the bulk with a very large energy 
density somewhere near the limiting parabola. It emits energy to the bulk, and after 
the initial decelerating period, it acquires an acceleration (presumably, the primordial 
"inflation"), which for appropriate choice of parameters and initial conditions can 
be of the order of 7 2 //5. It quickly cools down and starts decelerating, with energy 
density much smaller than the one it had initially. Assuming that the energy emission 
is very fast, it happens with the brane-bulk system out of equilibrium, and one may 
expect that the brane will supercool and reach a temperature much smaller than the 
roughly constant temperature of the bulk. As a consequence of the energy influx, 
the evolution of the brane will be attracted to the stable fixed point, analogous to 
the pi ^ discussed in the text. The corresponding ~ jpi -* will naturally be much 
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smaller than its early value, and may fit today's cosmic acceleration ||21|| . 

Clearly, the details of scenaria such as the above require further study and many 
open questions, such as the duration of accelerating periods, the creation of primor- 
dial fluctuations, and the compatibility with conventional cosmology at low energy 
densities, should be addressed. However, we believe that the cosmological evolution 
in the context we presented here has many novel features, that may provide answers 
to outstanding questions of modern cosmology. 

Some of the cosmological evolutions presented will have a (more general) holo- 
graphic dual that will be interesting to understand in more detail. This will provide 
a more controllable picture of the mechanisms advocated in this paper. 
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